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Abstract 

The power spectrum of local field potentials (LFPs) has been reported to scale as the inverse of 
CN ■ the frequency, but the origin of this "1// noise" is at present unclear. Macroscopic measurements 

^ . in cortical tissue demonstrated that electric conductivity (as well as permittivity) is frequency 

lO . dependent, while other measurements failed to evidence any dependence on frequency. In the 

present paper, we propose a model of the genesis of LFPs which accounts for the above data and 
contradictions. Starting from first principles (Maxwell equations), we introduce a macroscopic 
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O ' formalism in which macroscopic measurements are naturally incorporated, and also examine dif- 

QO ' 

^ ' ferent physical causes for the frequency dependence. We suggest that ionic diffusion primes over 

electric field eff"ects, and is responsible for the frequency dependence. This explains the contradic- 
• ^-j . tory observations, and also reproduces the 1// power spectral structure of LFPs, as well as more 

, complex frequency scaling. Finally, we suggest a measurement method to reveal the frequency 

I dependence of current propagation in biological tissue, and which could be used to directly test 

the predictions of the present formalism. 
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1 Introduction 



Macroscopic measurements of brain activity, such as the electroencephalogram (EEG), magnetoen- 
cephalogram or local field potentials (LFP), display approximately 1// frequency scaling in their 
power spectra (fTl l4l [30l [34l) . The origin of such 1// "noise" is at present unclear. 1// spectra can 
result from self-organized critical phenomena (|27]) . suggesting that neuronal activity may be work- 
ing according to such states @i. Alternatively, the 1/f scaling may be due to filtering properties of 
the currents through extracellular media dH). The latter hypothesis, however, was resting on indirect 
evidence, and still needs to be examined theoretically, which is one of the motivations of the present 
paper. 

A continuum model Q of LFPs incorporated the inhomogeneities of the extracellular medium into 
continuous spatial variations of conductivity (cr) and permittivity {e) parameters. This model repro- 
duced a form of low-pass frequency filtering in some conditions, while considering the extracellular 
medium as locally neutral with cr and e parameters independent of frequency. This model was not 
entirely correct, however, because macroscopic measurements in cortex revealed a frequency depen- 
dence of electrical parameters (fTSi) . We will show here that it is possible to keep the same model 
structure by including plausible causes for the frequency dependence. 

In a polarization model ^ of LFPs, the variations of conductivity and permittivity were considered 
by explicitly taking into account the presence of various cellular processes in the extracellular space 
around the current source. In particular, it was found that the phenomenon of surface polarization was 
fundamental to explain the frequency dependence of LFPs. The continuum model @ incorporated 
this effect phenomenologically through continuous variations of cr and e. In the polarization model, 
the extracellular medium is "reactive" in the sense that it reacts to the electric field by polarization 
effects. It is also locally non-neutral, which enables one to take into account the non-instantaneous 
character of polarization, which is at the origin of frequency dependence according to this model ([3]). 

In the present paper, we propose a "diffusion-polarization" model that synthesizes these previous ap- 
proaches and which takes into account both microscopic and macroscopic measurements. This model 
includes ionic diffusion, which we will show has a determinant influence on frequency filtering prop- 
erties. The model also includes electric polarization, which also influences the frequency-dependent 
electric properties of the tissue. We show that taking into account ionic diffusion and electric polar- 
ization allows us to quantitatively account for the macroscopic measurements of electric conductivity 
in cortical tissue according to the experiments of Gabriel et al. (fTSi) . 

However, recent measurements of Logothetis et al. (|^^ evidenced that the frequency dependence of 
cortical tissue was negligible, therefore in contradiction with the measurements of Gabriel et al. (fTSl) . 
We show here that the "diffusion-polarization" model can be consistent with both types of exper- 
iments, and thus may reconcile this contradiction. We will also examine whether this model can 
also explain the 1 // frequency scaling observed in LFP or EEG power spectra. Finally, we consider 
possible ways for experimental test of the predictions of this model. 

The final goal of this approach is to obtain a model of local field potentials which is consistent with 
both macroscopic measurements of conductivity and permittivity, and the microscopic features of the 
structure of the extracellular space around the current sources. 



2 Materials and Methods 



The numerical simulations presented in Section [33] describe the impedance of the extracellular medium 
based on the following equation 



1 r*oo 1 
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This equation gives the ciJ-frequency component of the impedance at point ri in extracellular space, 
in spherical symmetry (see ref. ^ and Eq.[lI]for details). 

To evaluate this equation, we use MATLAB, which computes the Riemann sum : 

Z(r co) = -y — ^- (2) 

where Ar' is the integration step (1 fim) and A'^ is determined for a "slice" of 1 mm. 



We also use the parametric model of Gabriel et al. (1151) to simulate the frequency dependence of 
electrical parameters cr and e of the extracellular fluid from gray matter (at a temperature of 37°C). 
This model is valid for frequencies included in the range of 10 Hz to 4 x 10^ Hz (fT6]) . According to 
this model, the absolute complex and macroscopic permittivity and conductivity (measured between 
10 and 10^*^ Hz) in cortical gray matter is given by the Cole-Cole parametric model (fTOl) : 



a> 



n=A 

y — — (3) 



.cr* Ae„ .cr 



where the sum runs over 4 polarization modes n, €„ = 8.85 x 10~^^ F/m is vacuum permittivity, foo = 
4.0 is the permittivity relative to / = 10^*^ Hz, cr = 0.02 S/m is the static conductivity at f = Hz 
according to the chosen parametric model. The parameters under the sum of Eq. [3] are given in 
Table 1. 



Table 1 here 



3 Results 



We start by outlining a macroscopic model with frequency-dependent electrical parameters (Sec- 
tion ISTT]), and we discuss the main physical causes for this frequency dependence (Section IX2l) . We 
then constrain the model to macroscopic measurements of electrical parameters, and provide numer- 
ical simulations to test the model and reproduce the experimental observations (Section IX^ . Finally, 
we propose a possible way to test the model experimentally (Section |34|) . 



3.1 Macroscopic model of local field potentials 

In this section, we derive the equations governing the time evolution of the extracellular potential. We 
follow a formalism similar to the one developed previously @, except that we reformulate the model 
macroscopically, to allow the electrical parameters (the conductivity cr and permittivity e) to depend 
on frequency, as demonstrated by macroscopic measurements (fT4l [T5l [T6l) . The physical causes of 
this macroscopic frequency dependence will be examined in Section IX2l 



3.1.1 General formalism 



We begin by deriving a general equation for the electrical potential when the electrical parameters are 
frequency dependent. We start from Maxwell equations, taking the 1 '''^ and the divergence of the 4^^'' 
Maxwell equation in a medium with constant magnetic permeability, giving: 

V ■ D = p-^''^'^ 

(4) 

where D, j and pJ''^'^ are respectively the electric displacement, current density and charge density in 
the medium surrounding the sources. 

Moreover, in a linear medium the equations linking the electric field ^ with electric displacement I), 
and with current density j, gives: 

e{x, t)^{x, t-T)dT (5) 

oo 

and 

a{x, t)E{x, t - T)dT . (6) 

oo 

The Fourier transforms of these equations are respectively Z)l, = e^^^l, and jX = o-^jEI), where we 
allow (T and e to depend on frequency. 

Given the limited precision of measurements, we can consider V xE »0 for frequencies smaller than 
1000 Hz. Thus, we can assume that ^ = - W such that the complex Fourier transform of Eqs.|4]can 
be written as: 

v-(6Uf)vy^) = -pi'"' 

Consequently, we have 

V ■ ((o-co + ioJ6c)VVaj) = V(cr^ + iaj6c) ■ + -((t^ + icoeJ^'^Va^ = (7) 



Compared to previous derivations (see Eq. 49 in ref. Q), this equation is a more general form in 
which the electrical parameters can be dependent on frequency. 



3.1.2 Macroscopic model 



In principle, it is sufficient to solve Eq. |7] in the extracellular medium to obtain the frequency de- 
pendence of LFPs. However, in practice, this equation cannot be solved because the structure of 
the medium is too complex to properly define the limit conditions. The associated values of elec- 
tric parameters must be specified for every point of space and for each frequency, which represents 
a considerable difficulty. One way to solve this problem is to consider a macroscopic or mean-field 
approach. This approach is justified here by the fact that the values measured experimentally are aver- 
aged values, which precision depends on the measurement technique. Because our goal is to simulate 
those measured values, we will use a macroscopic model, in which we take spatial averages of Eq.|7l 
and make a continuous approximation for the spatial variations of these average values. This type of 
approximation can be found in the classic theory of electromagnetism (see for example Chapters 9 
and 10 in Maxwell's treaty for electricity and magnetism 



To this end, we define macroscopic electric parameters, and cr^, as follows: 

=< eUx) >iy= fix,oj) 

and 

o'ffC^ =< (^oj(x) >\y= g{x,w) . 

where V is the volume over which the spatial average is taken. We assume that V is of the order of 
yum^, and is thus much smaller than the cortical volume, so that the mean values will be dependent of 
the position in cortex. 

Because the average values of electric parameters are statistically independent of the mean value of 
the electric field, we have: 

cr'^iT) < E >|, (f, t - T)dT + e^(r)— -^(f , t - T)dT , 

oo J —oo "I 

where the first term in the righthand represents the "dissipative" contribution, and the second term 
represents the "reactive" contribution (reaction from the medium). Here, all physical effects, such as 
diffusion, resistive and capacitive phenomena, are integrated into the frequency dependence of cr^ 
and e^. We will examine this frequency dependence more quantitatively in Section 1X21 

The complex Fourier transform of < f"'"' >\y {x, t) then becomes: 

< IT' >\v= (^ff + i^e^) < Ea. >w) = CT^ <E^>w , (8) 
where erf is the complex conductivity. We can also assume 

(rf = icoef (9) 

such that 

V- < t'"' >|.= V • (erf < E^%) = V • (icoef < eI >,,) = . (10) 

Because erf = (erf -i- zcuef ) and E^j >= -V < Vca >, the expressions above (Eqs. can also be 



written in the form: 



V-((cr^+iaje^)V<V^>ly) = 0. 



(11) 



We note that starting from the continuum model where only spatial variations were considered, 
and generalizing this model by including frequency-dependent electric parameters, gives the same 
mathematical form as the original model (compare with Eq 49 in ref. Q). This form invariance will 
allow us to introduce in Section [X2l the surface polarization phenomena by including an ad hoc fre- 
quency dependence in cr^ and 6^. The physical causes of this macroscopic frequency dependence 
is that the cortical medium is microscopically non-neutral (although the cortical tissue is macroscop- 
ically neutral). Such a local non-neutrality was already postulated in a previous model of surface 
polarization ([3]). This situation cannot be accounted by Eq.|7]if cr^ and 6^ are frequency independent 
(in which case p^, = when = 0). Thus, including the frequency dependence of these parameters 
enables the model to capture a much broader range of physical phenomena. 

Finally, a fundamental point is that the frequency dependences of the electrical parameters cr^ and 
cannot take arbitrary values, but are related to each-other by the Kramers-Kronig relations (fT3l 



where principal value integrals are used. These equations are valid for any linear medium (i.e., when 
Eqs.[5]and[6]are linear). These relations will turn out to be critical to relate the model to experiments, 
as we will see below. 

Note that, contrary to frequency dependence, the spatial dependences of crff and are independent 
of each other, because these dependences are related to the spatial distribution of elements within the 
extracellular medium. 

3.1.3 Simplified geometry for macroscopic parameters 

To obtain an expression for the extracellular potential, we need to solve Eq. dH which is possible 
analytically only if we consider a simplified geometry of the source and surrounding medium. The 
first simplification is to consider the source as monopolar. The choice of a monopolar source does 
not intrinsically reduce the validity of the results because multipolar configurations can be composed 
from the arrangement of a finite number of monopoles (35). In particular, if the physical nature of 
the extracellular medium determines a frequency dependence for a monopolar source, it will also 
do so for multipolar configurations. A second simplification will be to consider that the current 
source is spherical and that the potential is uniform on its surface. This simplification will enable 
us to calculate exact expressions for the extracellular potential and should not affect the results on 
frequency dependence. A third simplification is to consider the extracellular medium as isotropic. 
This assumption is certainly valid within a macroscopic approach, and justified by the fact that the 
neuropil of cerebral cortex is made of a quasi-random arrangement of cellular processes of very 
diverse size (81). This simplified geometry will allow us to determine how the physical nature of 
the extracellular medium can determine a frequency dependence of the LFPs, independently of other 




(12) 



and 




(13) 



factors (such as more realistic geometry, propagating potentials along dendrites, etc). 

Thus, considering a spherical source embedded in an isotropic medium with frequency-dependent 
electrical parameters, combining with Eq. [HI we have: 

(f <Voj>\y ^ 2 d< Vqj >iy ^ 1 d((Tca + iioeoj) d<Vaj>\y _^ ^^^^ 
dr^ r dr {(Toj + iojeoj) dr dr 

Integrating this equation gives the following relation between two points ri and rj in the extracellular 
space, 

2 d <V(jj>\y 2 d <Va) >\y 

^1 J (''l) [.cri^{ri) + ioje^^{ri)] = r2 (r2) [0-^(^2) + i(y6(^(r2)] . (15) 

^ dr ^ dr 

Assuming that the extracellular potential vanishes at large distances (< V^i >= 0), we find 

This equation is analogous to a similar expression derived previously (Eq. 25 in ref. dlj)), but more 
general. The two formalisms are related by: 

< loj >\v= CTz < ^OJ >\v= (O-^ + O < ^OJ >\v 

instead of jZ, = cr'^^oj (see Eq. 4 in ref. (2)). This difference is due to the fact that the conductivity 
here depends on frequency. 

In the following, we will use the simplified notations and V^j instead of < >|v,< E^oj >\v and 

< >\y, respectively. 

Using the relation V^j = Z(^Ico, the impedance Z^j is given by: 



1 

cr^inU+ioJT^ir')] 



where T(^(r') = and r is the distance between the center of the source and the position defined by 



3.2 Physical causes for frequency-dependent electrical parameters 

In the following, we successively consider two different cases of extracellular medium: first, "non- 
reactive" media, in which the current passively flows into the medium; second, "reactive" media, in 
which some properties (such as charge distribution) may change following current flow. For each 
medium, we will consider two types of physical phenomena, the current produced by the electric 
field, and the current produced by ionic diffusion, as schematized in Fig.[T]). 



3.2.1 Non-reactive media with electric fields (Model N) 



"Non reactive" media (^^%- << 1 , crf^ = cr^ and 6^ = e^) are equivalent to "resistive" media, in which 
the resistance (or equivalently, the conductivity) does not change following the flow of current. The 
simplest type of such configuration consists of a resistive medium (such as a homogeneous conduc- 
tive fluid) in which current sources solely interact via their electric field. Applying Eq[l7] to this 
configuration is equivalent to model the extracellular potential by Coulomb's law: 

where Va,(r) is the extracellular potential at a position defined by r in extracellular space, r is the 
absolute distance between r and the center of the current source. Here, the conductivity (cr^(r) = cr^) 
is independent of space and frequency, and thus, this model is not compatible with macroscopic 
measurements of frequency dependence (fT4l[T5l[T6l) . It is, however, the most frequently used model 
to calculate extracellular field potentials (|3TI) . This model will be referred to as "Model N" in the 
following. 



3.2.2 Non-reactive media with ionic diffusion (Model D) 

Because current sources are ionic currents, there is flow of ions inside or outside of the membrane, 
and another physical phenomena underlying current flow is ionic difl'usion. Let us consider a resistive 
medium such as a homogeneous extracellular conductive fluid with electric parameters 

and in which the ionic difl'usion coefficient is D. When the extracellular current is exclusively due 
to ionic diffusion, the current density depends on frequency as ^/cij (see Appendix |B]). A resistive 
medium behaves as if it had a resistivity equal to a(l + -^), where b is complex. The parameter 
a is the resistivity for very high frequencies, and reflects the fact that the effect of ionic diffusion 
becomes negligible compared to calorific dissipation (Ohm's law) for very high frequencies. When 
ionic diffusion is dominant compared to electric field effects, the real part of b is much larger than a. 

The frequency dependence of conductivity will be given by: 

= ^ (19) 
yo) +k 

Applying Eq. [l7] to this configuration gives the following expression for the electric potential as a 
function of distance: 

^'^'^^ = ~A — M ■ ~ = — ■ — m ■ ~ 



This expression shows that, in a non-reactive medium, when the extracellular current is dominated 
by ionic diffusion (compared to that directly produced by the electric field), then the conductivity 
will be frequency dependent and will scale as ^f(ji>. This model will be referred to as "Model D" in 



the following. Note that, if the electric field primes over ionic diffusion, then we have the situation 
described by Model N above. 



3.2.3 Reactive media with electric fields (Model P) 



In reality, extracellular media contain different charge densities, for example due to the fact that cells 
have a non-zero membrane potential by maintaining differences of ionic concentrations between the 
inside and outside of the cell. Such charge densities will necessarily be influenced by the electric field 
or by ionic diffusion. As above, we first consider the case with only electric field effects and will 
consider next the influence of diffusion and the two phenomena taken together. 

Electric polarization is a prominent type of "reaction" of the extracellular medium to the electric field. 
In particular, the ionic charges accumulated over the surface of cells will migrate and polarize the cell 
under the action of the electric field. It was shown previously in a theoretical study that this "surface 
polarization" phenomena can have important effects on the propagation of local field potentials (3). If 
a charged membrane is placed inside an electric field Eq , there is production of a secondary electric 
field ^f, given by (see Eq. 31 in ref. 

€ = 7—^- (21) 

1 + l(OTM 



This expression is the frequency-domain representation of the effect of the inertia of charge movement 
associated with surface polarization, reflecting the fact that the polarization does not occur instanta- 
neously but requires a certain time to setup. This frequency dependence of the secondary electric 
field was derived in ref. (Sj) for a situation where the current was exclusively produced by electric 
field. The parameter tm is the characteristic time for charge movement (Maxwell- Wagner time) and 
equals j ^memb ^ where e'"^'^'^ and cr"'^'"'^ are respectively the absolute (tangential) permittivity 
and conductivity of the membrane surface, respectively, and are in general very different from the 
permittivity and the conductivity of the extracellular fluid. 

Let us now determine for zero-frequency the amplitude of the secondary field produced between 
two cells embedded in a given electric field. First, we assume that it is always possible to trace a con- 
tinuous path which links two arbitrary points in the extracellular fluid (see Fig. |2^). Consequently, 
the domain defined by extracellular fluid is said to be linearly connex. In this case, the electric po- 
tential arising from a current source is necessarily continuous in the extracellular fluid. Second, in a 
first approximation, we can consider that the cellular processes surrounding sources are arranged ran- 
domly (by opposition to being regularly structured) and their distribution is therefore approximately 
isotropic. Consequently, the field produced by a given source in such a medium will also be ap- 
proximately isotropic. Also consequent to this quasi-random arrangement, the equipotential surfaces 
around a spherical source will necessarily cut the cellular processes around the source (Fig. [2^). 

Suppose that at time ? = 0, an excess of charge appears at a given point in extracellular space, then 
a static electric field is immediately produced. At this time, currents begins to flow in extracellular 
fluid, as well as inside the different cellular processes surrounding the source. These cells begin to 
polarize, with tmw as the characteristic polarization time. Asymptotically, the system will reach an 
equilibrium where the polarization will neutralize the electric field, such that there is no electric field 



inside the cells (zero current). Now suppose that, in the asymptotic regime, there would still be a 
current flowing in between cells (in the extracellular fluid), then we have two possibilities. First, the 
equipotential surfaces are discontinuous, or they cut the membrane surfaces (as illustrated in Fig. [2^). 
The first possibility is impossible because it would imply an infinite electric field. The second pos- 
sibility is also impossible, because cells are isopotential due to polarization. Therefore, we can say 
that, asymptotically, there is no current flowing in extracellular fluid at / = 0, and necessarily this is 
equivalent to a dielectric medium. In other words, a passive inhomogeneous medium with randomly 
distributed passive cells is a perfect dielectric at zero frequency . In this case, the conductivity must 
tend to when frequency tends to 0. Thus, in the following, we assume that = -E^ where E^ is 
the field produced by the source. 



Figure |2] here 

It follows that the expression for the current density in extracellular space as a function of the electric 
field is given by: 

L = crfi'^ = CT"^ ■ (1+/—) • (r + £^^) = (r'«. (1+/—) • • 



In addition, for cerebral cortical tissue, we have 1 + ~ 1 for frequencies > 10 Hz and < 1000 Hz 
(see ref. (fT5l) . Thus, an excellent approximation of the conductivity can be written as: 



M^^..^£^ (22) 

1 + itOTM 



Applying Eq.fTT] gives: 



ya>{r) = TT-— = -— — -■— (23) 

AncTy r 1 + icotm Ancr'" r 



This model describes the eff"ect of polarization in reaction to the source electric field, and will be 
referred to as "Model P' in the following. 



3.2.4 Reactive media with electric field and ionic diff"usion (Model DP) 



Figure [H here 

The propagation of current in the medium is dominated by ionic diff"usion currents or by currents 
produced by the electric field, according to the values of k and ki with respect to ^/a7. The values 
of k and ki are respectively inversely proportional to the square root of the global ionic difl^usion 
coefficient in the extracellular fluid, and of membrane surface (see Appendix |B]). 



We apply the reasoning based on the connex topology of the cortical medium (see above) to deduce 



the order of magnitude of the induced field for zero frequency 



1 + 1 yo) T 



where 



^memb 

T = (yfo7 + kl)TM = Vtt» 



0-memb 



Because the "tangential" conductivity on membrane surface is given by 

^Memb . r~~ 
memb ^ V"^ 



yfto +ki 

when the current is dominated by either electric field or ionic diffusion (see Eq. \T% . 

It follows that the expression for the current density in extracellular space as a function of the electric 
field is given by 

yo) + k cr cu yoj + k l+i yo) t 

because 1 + i—M- cortical tissue for frequencies > 10 Hz and < 1000 Hz (see ref. (fTSl) ). 
Thus, we have the following expression for the complex conductivity of the extracellular medium: 

ya> + k 1 +1 ya> r 

where r = ( yfoj + ki)TM- 

Thus, we have obtained a unique expression (Eq. [251) for the apparent conductivity in extracellu- 
lar space outside of the source, and its frequency dependence due to differential Ohm's law, elec- 
tric polarization phenomena and ionic diffusion. These phenomena are responsible for an apparent 
frequency-dependence of the electric parameters, which will be compared to the frequency depen- 
dence observed in macroscopic measurements of conductivity (Section [13]) . 

Finally, Eqs. [IT] and [25] imply that the macroscopic impedance of a homogeneous spherical shell of 
width R2-R1 is given by: 

7=1 r ^ _ i?2-i?l 1 .26) 



In the following, this model will be referred as the "diffusion-polarization" model, or "DP" model, 
and we will use the above expressions (Eqs. [25] and [26]) to simulate experimental measurements. 



3.3 Numerical simulation of macroscopic measurements 



3.3.1 Experiments of Gabriel et al. (1996) 



We first consider the experiments of Gabriel et al. (fT?l [T5l [T6l) who measured the frequency depen- 
dence of electrical parameters for a large number of biological tissues. In these experiments, the 
biological tissue was placed in between two capacitor plates, which were used to measure the ca- 
pacitance and leak current using the relation /^^ = YVoj, imposing the same current amplitude at all 
frequencies. Because the admittance value is proportional to crf^ + i(x>e^ , measuring the admittance 
provides direct information about cr^ and . 

To stay coherent with the formalism developed above, we will assume that the capacitor has a spher- 
ical geometry. The exact geometry of the capacitor should in principle have no influence on the 
frequency dependence of the admittance, because the geometry will only aff"ect the proportionality 
constant between cr^ and 7^. In the case of a spherical capacitor, by applying Eql26l we obtain: 

K K2—K1 K2—K1 



We also take into account the fact that the resistive part is always greater than the reactive part for low 
frequencies (< 1000 Hz), which is expressed by 

^^ijj I^^OJ «^ ■ 

This relation can be verified for example from the Gabriel et al. measurement (TS), where it is true for 
the whole frequency band investigated experimentally (between 10 and 10^*^ Hz). 



The real part of cr^ = cTy then takes the form 



cr,. 



.M 



cr 



M 



OJT 



y/co + k (x>T^ + 1 
where r = ( + ^1 )tm 

By substituting this value of t, the inverse of the conductivity (the resistivity) is given by: 



(28) 



1 Ilk 

1 + ■ 
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CO + 2k\ -^Jco + fcj ^ 



2 
1 

(29) 



Finally, to reproduce Gabriel et al. experiments, we assume ki » ^[a>. By developing in series the 
last term (in parenthesis) of Eq.[29l we obtain: 



\ — K\ K2 K3 Ki K2 K3 

— ~Ko + —r^ + — + = Ko + —r^ + — + ■ 



cr. 



CO^I^ CO CO 



,3/2 



fill f pi2 



(30) 



Figure [3] here 



Eq. [30] corresponds to the conductivity cr^, as measured in the experimental conditions of Gabriel 
et al. experiments (the permittivity is obtained by applying Kramers-Kronig relations). Figure 
shows that this expression for the conductivity can explains the measurements in the frequency range 
of 10 to 1000 Hz, which are relevant for LFPs. To obtain this agreement, we had to assume in Eq.[25]a 
relatively low Maxwell- Wagner time of the order of 0.15 s (fc = 1 /(Itttm) between 1 Hz and 10 Hz), 
ki > ^[oj > k (for frequencies smaller than 100 Hz). 

This value of Maxwell-Wagner time is necessary to explain Gabriel's experiments, and may seem 
very large at first sight. However, the Maxwell- Wagner time is not limited by physical constraints, 
because we have by definition tmw = In principle, the value of cr^ can be very small, approaching 
zero, while e^j can take very large values. For example, taking Gabriel et al. measurements in aqueous 
solutions of NaCl and in gray matter (fT5l) . gives values of tmw comprised between 1 ms and 100 ms 
for frequencies around 10 Hz. 

Thus, the model predicts that in the Gabriel et al. experiments, the transformation of electric current 
carried by electrons to ionic current in the biological medium necessarily implies an accumulation 
of ions at the plates of the capacitor. This ion accumulation will in general depend on frequency, 
because the conductivity and permittivity of the biological medium are frequency dependent. This 
will create a concentration gradient across the biological medium, which will cause a ionic diff"usion 
current opposite to the electric current. This ionic current will allow a greater resulting current be- 
cause surface polarization is opposite to the electric field. Figure O shows that such conditions give 
frequency-dependent macroscopic parameters consistent with Gabriel et al. measurements. 

The parameter choices to obtain this agreement can be justified qualitatively because the ionic diffu- 
sion constant on cellular surfaces is probably much smaller than in the extracellular fluid, such that 
k\ » k. This implies the existence of a frequency band Bf for which is negligible with respect 
to k\ , but not with respect to k because these constants are inversely proportional to the square root of 
their respective diffusion coefficients. Thus, the approximation that we suggest here is that this band 
Bf finishes around 100 Hz in the Gabriel et al. experimental conditions. It is important to note that 
this parameter choice is entirely dependent on the ratio between ionic diffusion current and the current 
produced by the electric field, and thus will be depend on the particular experimental conditions. 

It is interesting to note that the present model and the phenomenological Cole-Cole model (|TOl) pre- 
dict different behaviors of the conductivity for low frequencies (<10 Hz). In the present model, the 
conductivity tends to zero when frequency tends to zero, while in the Cole-Cole extrapolation, it tends 
to a constant value (fT4l) . The main difference between these models is that the Cole-Cole model is 
phenomenological and has never been deduced from physical principles for low frequencies, unlike 
the present model which is entirely deduced from well-defined physical phenomena. 

3.3.2 Logothetis et al. (2007) measurements 

We next consider the experiments of Logothetis et al. (|^^ . which reported a resistive medium, in 
contrast with Gabriel et al. experiments. In these experiments, 4 electrodes were aligned and spaced 
by 3 mm in monkey cortex. The first and last electrodes were used to inject current, while the two 
intermediate electrodes were used to measure the extracellular voltage. The voltage was measured at 
different frequencies and current intensities. 



One important point in this experimental setup is that the intensity of the current was such that the 
vokage at the extreme (injecting) electrodes saturates. One of the consequences of this saturation was 
to limit ionic diffusion effects, as discussed by the authors (1231) . This voltage saturation will diminish 
the concentration difference near the source (we would have an amplification if this was not the case). 
It follows that, in the Logothetis et al. experiments, the ratio between diffusion current and electric 
field current is very small. Thus, in this case we use values of parameters ki very different from those 
assumed above to reproduce Gabriel et al. experiments, in particular ki « y/oj. As we will see in the 
Discussion section, this situation may be different from the physiological conditions. 

Nevertheless, the large distance between electrodes suggests that the relation between current and 
voltage is linear because the current density is roughly proportional to the inverse of squared distance 
to the source. Consequently, we can suppose that in the Logothetis et al. experiments, the ionic 
gradient is negligible, which prevents ionic diffusion currents. Thus, in this experiment, most - if not 
all - of the extracellular current is due to electric-field effects. 



In this situation, the conductivity (Eq.[25l) becomes: 



which is similar to the "Model P" above. 



Moreover, taking the same Maxwell-Wagner time tm as above for Gabriel et al. experiments (which 
corresponds to a cutoff frequency of 1 Hz), we have for frequencies above 10 Hz: 



(iOTMr _ 



similar to a resistive medium. 

Thus, in the Logothetis et al. experiments, the saturation phenomenon entrains current propagation 
in the biological medium as if the medium was quasi-resistive for frequencies larger than 10 Hz. 
This constitutes a possible explanation of the contrasting results in Logothetis et al. and Gabriel et al. 
measurements. 



3.3.3 Frequency dependence of the power spectral density of extracellular potentials 

The third type of experimental observation is the fact that the power spectral density (PSD) of LFPs 
or EEG signals displays 1/f frequency scaling ([H |4l [301 1341) . To examine if this 1/f scaling can be 
accounted for by the present formalism, we consider a spherical current source embedded in a contin- 
uous macroscopic medium. We also assume that the PSD of the current source is a Lorentzian, which 
could derive for example from randomly occurring exponentially decaying postsynaptic currents (H)) 
(see Fig. HI). 

To simulate this situation we used the "diffusion-polarization" model with ionic diffusion and electric 
field effects in a reactive medium. We have estimated above that surface polarization phenomena have 
a cutoff frequency of the order of 1 Hz, and will not play a role above that frequency. So, if we focus 
on the PSD of extracellular potentials in the frequency range larger than 1 Hz, we can consider only 



the effect of ionic diffusion (in agreement with the Gabriel et al. experiments - see above). 
Thus, we can approximate the conductivity as (see Eq.[25l): 

(T^ = ayfZj (32) 

where a is a constant. 

It follows that the extracellular voltage around a spherical current source is given by (see Eg. [171): 

Vir,a)) = — = — — = — (33) 

4na yio r yco r yo) 

where R is the radius of the source. 

Figure HI here 

In other words, we can say that the extracellular potential is given by the current source Icu convolved 
with a filter in 1/ y/co, which is essentially due to ionic diffusion (Warburg impedance; see refs. ([T2l 
[I9l|20l)). A white noise current source will thus result in a PSD scaling as 1/f, and can explain the 
experimental observations, as shown in Fig. |4l Experimentally recorded LFPs in cat parietal cortex 
display LFPs with frequency scaling as 1/f for low frequencies, and 1/f^ for high frequencies (Fig.|3]\- 
B). Following the same procedure as in ref. (4), we reconstructed the synaptic current source from 
experimentally recorded spike trains (Fig.|4lC-D). The PSD of the current source scales as a Lorentzian 
(Fig. 13^) as expected from the exponential nature of synaptic currents. Calculating the LFP around 
the source and taking into account ionic diffusion, gives a PSD with two frequency bands, scaling in 
1/f for low frequencies, and 1/f^ for high frequencies (Fig.|3^). This is the frequency scaling observed 
experimentally for LFPs in awake cat cortex dlj). We conclude that ionic diffusion is a plausible 
physical cause of the 1 // structure of LFPs for low frequencies. 

Two important points must be noted. First, the diffusion-plolarization model does not automatically 
predict l/f scaling at low frequencies, but rather implements an 1// filter, which may result in fre- 
quency scaling with larger slopes. Second, the same experimental situation may result in different 
frequency scaling, and this is also consistent with the diffusion-polarization model. These two points 
are illustrated in Fig.[5]which shows a similar analysis as Fig.|4]but during slow-wave sleep in the same 
experiment. The LFP is dominated by slow-wave activity (Fig. [IK), and the different units display 
firing patterns characterized by concerted "pauses" (gray lines in Fig. [5^), characteristic of slow-wave 
sleep and which are also visible in the reconstructed synaptic current (Fig.[5]C!). The PSD shows a sim- 
ilar scaling as 1//^ as for wakefulness, but the scaling at low frequencies is different (slope around 
-2 at low frequencies; see Fig. [SjD). The PSD reconstructed using the diffusion-polarization model 
displays similar features (compare with Fig. [5^). This analysis shows that the diffusion-polarization 
model qualitatively accounts for different regions of frequency scaling found experimentally in dif- 
ferent frequency bands and network states. 



Figure [5] here 



3.4 Measurement of frequency dependence 



In this final section, we examine a possible way to test the model experimentally. The main prediction 
of the model is that, in natural conditions, the extracellular current perpendicular to the source is 
dominated by ionic diffusion. The experiments realized so far (fT4l USl IT6l 123]) used macroscopic 
currents that did not necessarily respect the correct current flow in the tissue. 

We suggest to create more naturalistic current sources by creating ionic currents using a micropipette 
placed in the extracellular medium. By using periodic current injection during very short time A? 
compared to the period (small duty cycle), we can measure using the same electrode the extracellular 
voltage Vpr (using a fixed reference far away from the source). If the period of the source is shorter 
than the relaxation time of the system, the voltage will integrate, which is due to charge accumulation. 

Because Vpr is directly proportional to the amount of charge emitted as a function of time during At 
(capacitive effect of the extracellular medium), the time variation of Vpr is directly proportional to the 
ionic diffusion current. In such conditions, if the extracellular medium is purely resistive, as predicted 
by Logothetis et al. experiments, the relaxation time should be very small, of the order of 10~^^ s 
@, which would prevent any integration phenomena and charge accumulation for frequencies below 
10^^ Hz. If the medium has a slower relaxation, due to polarization and ionic diffusion, then we should 
observe voltage integration and charge accumulation for physiological frequencies (< 1000 Hz). 

To illustrate the difference between these two situations, we consider the simplest case of a non- 
reactive medium (as in Model D above), in which the current can be produced by ionic diffusion or by 
the electric field, or by both. To calculate the time variations of ionic concentration and extracellular 
voltage, we consider the current density: 

j = DVe[c]+(rE (34) 
According to the differential law for charge conservation and Poisson law, we have: 

V/+|=DVV+^P+|=0 (35) 

ot e ot 

When ionic diffusion is negligible compared to Ohm's law, we have: 

dp cr 

f + -p«0 

ot e 

It follows that the charge density is given by: 

p=Poexp^-^fj (36) 

On the other hand, if ionic diffusion is the primary cause of current propagation in the extracellular 
medium, then the relaxation time should be much larger and thus, integration should be observed. 
When ionic diffusion is dominant, we have: 

f + DV2p.O 
at 



instead of Eq.[35i 



The general solution is: 

p = —== p(r,0)e~^dr (37) 
y/lDt Jo 

The difference between the expressions above (Eqs.[36land[T7]) shows that the time variation of charge 
density is different according to which current dominates, electric field current or ionic diffusion cur- 
rent. The same applies to the electric potential between the electrode and a given reference, because 
this potential is linked to charge density through Poisson's law. Therefore, this experiment would be 
crucial to clearly show which of the two current primes for currents perpendicular to the source (this 
would not apply to longitudinal currents, like axial currents in dendrites). In the hypothetical case of 
dominant ionic diffusion, the cortex would be similar to a Warburg impedance and one can estimate 
the macroscopic diffusion coefficient using Eq. [J7i 

Thus, using a micropipette injecting periodic current pulses, it should be possible to test the capacity 
of the medium to create charge accumulation for physiological frequencies. If this is the case, this 
would constitute evidence that ionic currents are non-negligible in the physiological situation. 



4 Discussion 



In the present paper, we have proposed a framework to model local field potentials, and which syn- 
thesizes previous measurements and models. This framework integrates microscopic measurements 
of electric parameters (conductivity cr and permittivity e) of extracellular fluids, with macroscopic 
measurements of those parameters (cr^, e^) in cortical tissue (fTSl E^ . It also integrates previous 
models of LFPs, such as the continuum model ©, which was based on a continuum hypothesis of 
electric parameters variations in extracellular space, or the polarization model ([3]), which explicitly 
considered different media (fluid and membranes) and their polarization by the current sources. The 
present model is more general and also integrates ionic diffusion, which is predicted as a major deter- 
minant of the frequency dependence of LFPs. This "diffusion-polarization" model also accounts for 
observations of 1// frequency scaling of LFP power spectra, which is due here to ionic diffusion, and 
is therefore predicted to be a consequences of the genesis of the LFP signal, rather than being solely 
due to neuronal activity (see ref. (4)). Finally, this work suggests that ephatic interactions between 
neurons can occur not only through electric fields but that ionic diffusion should also be considered 
in such interactions. 

As discussed in Section 13.1.31 the present model rests on several approximations, which were nec- 
essary to obtain the analytic expressions used here. These approximations were that current sources 
were considered as monopolar entities (longitudinal currents such as axial currents in dendrites were 
not taken into account), the current source was spherical and the extracellular medium was considered 
isotropic. Because multipolar effects can be reconstructed from the superposition of monopoles (|35l) . 
the monopolar configuration should not affect the results on frequency dependence, as long as the 
extracellular current perpendicular to the source is considered. Similarly, the exact geometry of the 
current source should have no influence on the frequency dependence far away from the sources. 
However, in the immediate vicinity of the sources, the geometric nature and the synchrony of synaptic 



currents can have influences on the power spectrum (|33l) . Another assumption is that the extracellular 
medium is isotropic, which was justified within the macroscopic framework followed here. These 
factors, however, will influence the exact shape of the LFPs. More quantitative models including a 
more sophisticated geometry of current sources and the presence of membrane excitability and action 
potentials should be considered (e.g., see refs. (fT8ll33]) ). 

The main prediction of the present model is that ionic diff"usion is an essential physical cause for the 
frequency dependence of LFPs. We have shown that the presence of ionic diffusion allows the model 
to account quantitatively for the macroscopic measurements of the frequency dependence of electric 
parameters in cortical tissue flST) . Ionic diffusion is responsible for a frequency dependence of the 
impedance as 1 / ^fo) for low frequencies (< 1000 Hz), which directly accounts for the observed 1 // 
frequency scaling of LP? and EEG power spectra during wakefulness ([Ill4l[30l|34l) (see Fig.|4l). Note 
that the EEG is more complex because it depends on the diffusion of electric signals across fluids, 
dura matter, skull, muscles and skin. However, this "filtering" is of low-pass type, and may not affect 
the low-frequency band, so there is a possibility that the 1/f scaling of EEG and LFPs have a common 
origin. The present model is consistent with the view that this apparent "1// noise" in brain signals 
is not generated by self-organized features of brain activity, but is rather a consequence of the genesis 
of the signal and its propagation through extracellular space Q. 

It is important to note that the fact that ionic diffusion may be responsible for 1 // frequency scal- 
ing of LFPs is not inconsistent with other factors which may also influence frequency scaling. For 
example, the statistics of network activity - and more generally network state - can affect frequency 
scaling. This is apparent when comparing awake and slow-wave sleep LFP recordings in the same 
experiment, showing that the 1 // scaling is only seen in wakefulness but 1 //^ scaling is rather seen 
during sleep dH) (see Fig. [5]). In agreement with this, recent results indicate that the correlation struc- 
ture of synaptic activity may influence frequency scaling at the level of the membrane potential, and 
that correlated network states scale with larger (more negative) exponents (|25|) . 

We also investigated ways to explain the measurements of Logothetis et al. (23), who reported that 
the extracellular medium was resistive and therefore did not display frequency dependence, in con- 
tradiction with Gabriel et al. (fTSi) measurements. We summarize and discuss our conclusions below. 

In Gabriel et al. experiments (IT5] ). one measures permittivity and conductivity in the medium in be- 
tween two metal plates. This forms a capacitor, which (macroscopic) complex impedance is mea- 
sured. This measure actually consists of two independent measurements, the real and imaginary part 
of the impedance. These values are used to deduce the macroscopic permittivity and macroscopic 
conductivity of the medium. However, at the interface between the medium and the metal plates, the 
flow of electrons in the metal corresponds to a flow of charges in the tissue, and a variety of phenom- 
ena can occur, which can interfere the measurement. The accumulation of charges that occurs at the 
interface between the electrode and the extracellular fluid implies a polarization impedance, which 
depends on the interaction between ions and the metal plate. Because this accumulation of charge 
implies a variation of concentration, the flow of ions may involve an important component of ionic 
diffusion. 

In Logothetis et al. experiments, a system of 4 electrodes is used, the two extreme electrodes inject 
current in the medium, while the two electrodes in the middle are exclusively used to measure the volt- 
age. This system is supposed to be more accurate than Gabriel's, because the electrodes that measure 
voltage are not subject to charge accumulation. However, the drawback of this method are nonlin- 



ear effects. The magnitude of the injected current is such that the voltage at the extreme electrodes 
saturates. This voltage saturation also implies saturation of concentration (capacitive effect between 
electrodes), which limits ionic diffusion currents. Thus, the ratio between ionic diffusion currents and 
the currents due to the electric field is greatly diminished relative to Gabriel et al. experiments. 

We think that natural current sources are closer to Gabriel et al. situation for several reasons. First, 
the magnitude of the currents produced by biological sources is far too low for saturation effects. 
Second, the flow of charges across ion channels will produce perturbations of ionic concentration, 
which will be re-equilibrated by diffusion. The effects may not be as strong as the perturbations of 
concentrations induced by Gabriel et al. type of experiments, but ionic diffusion should play a role in 
both cases. This is precisely one of the aspects that should be evaluated in further experiments. 

The Logothetis et al. experiments were done using a 4-electrode setup which neutralizes the influ- 
ence of electrode impedance on voltage measurements (fTTl l26l) . This system was used to perform 
high-precision impedance measurements, also avoiding ionic diffusion effects (1231) . Indeed, these ex- 
perimental conditions, and the apparent resistive medium, could be reproduced by the present model 
if ionic diffusion was neglected. The present model therefore formulates the strong prediction that 
ionic diffusion is important, and that any measurement technique should allow ionic diffusion to re- 
veal the correct frequency-dependent properties of impedance and electric parameters in biological 
tissue. 

The critical question that remains to be solved is whether, in physiological conditions, ionic diffusion 
plays a role as important as suggested here. We propose a simple method to test this hypothesis. 
The frequency dependence could be evaluated by using an extracellular electrode injecting current in 
conditions as close as possible to physiological conditions (a micropipette would be appropriate). By 
measuring the integration of the extracellular voltage following periodic current injection, one could 
estimate the "relaxation time" of the medium with respect to charge accumulation. If this relaxation 
time occurs at time scales relevant to neuronal currents (milliseconds) rather than the fast relaxation 
predicted by a purely resistive medium (picoseconds), then ionic diffusion will necessarily occur in 
physiological conditions, which would provide evidence in favor of the present mechanism. 



Appendices 

A Estimation of ionic diffusion current vs. electric diffusion in 
sea water and in cortex. 

In this appendix, we evaluate the ratio between ionic diffusion currents and electric field currents 
in the extracellular space directly adjacent to the source. This ratio measures if the ionic diffusion 
current perpendicular to the membrane is greater than the electric field current. We will design this 
ratio by the term r/g. 



We have in general 



iTotal = eDV{C]+(TVV , 



(38) 



where the first term in right hand is the electric current density produced by ionic diffusion, while the 
second term is that produced by differential Ohm's law. 



For a displacement Ar = 10 nm in the direction across the membrane (from inside to outside), we have 
approximatively : 

hotal ■ dr ^ jTotal ' Ar = eA[C]^r=lQnm + O-ca/^V^r=l0mn (39) 

Suppose that we have a spherical cell of 10 fim radius, embedded in sea water and at resting potential. 
The resting membrane potential is a dynamic equilibrium between inflow and outflow of charges, in 
which these two fluxes are equal on (temporal) average. Fluctuations of current around this average in 
the extracellular medium around the membrane have all characteristics of thermal noise (l32l) because 
the shot noise (see refs. @ and (|38] )) is zero when the current is zero on average, such that the net 
charge on the external side of the membrane varies around a mean value with the same characteristics 
as white noise (thermal noise). These fluctuations will therefore be present also at the level of the 
membrane potential. In this appendix, we evaluate the order of magnitude of the electric current 
caused by ionic diffusion, relative to the electric field for this situation of dynamic equilibrium. 

First, the ratio between the membrane voltage noise and the variation of total charge concentration is 
given by 

AQtot — C AVfiiembrane — ^1^^ membrane — 1-25 X 10 AVuigmhi-ane > (40) 

because the cell's capacitance is given by C = AnR^Cm = O.OAnR^ (Cm - 10"^ FjnP-) 
Second, mass conservation imposes: 

AQtot = e*Veff*/S\C]tot , (41) 

where e = 1 .69 x 10" C and Veff is the volume of the spherical shell containing the charges. Because 
the charges are not uniformly distributed inside the cell, but rather distributed within a thin spherical 
shell adjacent to the membrane, because the electric field developed across the membrane is very 
intense (of the order of V/m = 10^ V/m), Thus, the width of the shell is of the order of 

dR - IQ-'^R < 1 nm, where the volume of the spherical shell is approximately equal to AnR^dR. In 
this case, we have for monovalent ions (|z| = 1) 

^Qto, = k2A[C]tot - 2.2 X 10-38a[C],o, = 2.2 X lO-^^A[C]^r=10nm (42) 

if we assume that the variation of concentration on the adjacent border of the exterior surface of the 
cell is over a width of 10 nm. 

In this case, we have 

A[C]a,= io„„, ^ |i ^ ^q27 ^/,^3^ (43) 
^ 'membrane i^l 

Third, the potential difference between the cell surface and 10 nm away of it, is given by 

AVa r=\Qnm ~ membrane „ — 10 AVmgmbrane 

(44) 

R + r 



Consequently, the ratio between ionic diffusion current and electric diffusion current caused by ther- 



mal noise in sea water obeys 



fie ~ ^sea \-\7 ~ ^seai, ~ ^sea ' ^ ' 



where the diffusion constant of or Na'^ in sea water is of the order of 10 ^ This implies that 
the ratio r is much larger than 1 for frequencies <1000 Hz because cr^^" of sea water is necessarily 
<2S/m. 



Because tortuosity is given by A = ^ D sea /D cortex, and is comprised between 1.6 and 2.2 (for small 
and large molecules, respectively) in cerebral cortex (|27l[28l[36l) . the macroscopic diffusion constant 
in cortex is certainly larger than Dsea/^0. Thus, we have 



_10 



C"" > = (46) 



where crf^''^ ^_ ^OAS/m (see ref. (TF)). 

This evaluation shows that the phenomenon of ionic diffusion is essential to determine the current 
field in the cortex. 

Finally, we note that we did not need to evaluate the absolute magnitude of AV in our evaluation. This 
evaluation is valid for a physical situation where we have a permanent white noise over a distance of 
10 nm, independently of the intensity of this noise (which in practice will depend on many factors, 
such as the size of the cell, the number of ion channels, etc). 



B Frequency scaling of ionic diffusion 

In this appendix, we calculate the frequency dependence of ionic diffusion current outside of a spher- 
ical current source. We consider a constant variation of ionic concentration, AX,, on the surface of the 
source, and a null variation at an infinite distance (Warburg conditions). 

The diffusion equation for a given ionic species is 

— -^=AV'AX,, (47) 
ot 

where AX, is the perturbation of concentration X,- of ion / around the steady-state value, and D,- is the 
associated diffusion coefficient. This diffusion coefficient depends of the ionic species considered and 
the structure of the medium. 

Because the geometry of the problem and the limit conditions respect spherical symmetry, we use 
spherical coordinates. In this coordinate system, we have 

(48) 

ot dr^ r or 



because AX, does not depend on 6 and of O (spheric symmetry). 



The Fourier transform of AX, with respect to time gives: 



" + = ^ + = —AXi (49) 



dr^ r dr dr^ r dr Di 
which general solution is given by 



AX, = AM + (50) 



For a variation of concentration at the source border which is independent of frequency and which 
satisfies the Warburg hypothesis (the variation of concentration tends to zero at an infinite distance (fT^ 
[37])), we have: 

AX,-Jr) = AX,-Ji?) • (51) 

where r is the distance between the center of the source and R is the radius of the source. 

Thus, the electric current density produced by ionic diffusion is given by: 

ji(r) = ZeDi—^f = -ZeDi ( - + J— ) AX,Jr) r (52) 
or r y Di 

where Ze is the charge of ions /. 

Because we can consider that the source and extracellular medium form a spherical capacitor, the 
voltage difference between the surface of the source and infinite distance is given by ZeCAXi^(R) 
where C is the capacitance value. Thus, the electric impedance of the medium is given by: 

= (53) 



For a source of radius R = 10 /um and a macroscopic ionic diffusion coefficient of the order of 
10~^^ m^/s, and for frequencies > 1 Hz, we can approximate the impedance by: 

Zoj « (54) 



The same expression for the impedance is also obtained in cylindric coordinates or planar Cartesian 
coordinates (not shown). 

Note that if several ionic species are present, then the superposition principle applies (Fick equations 
are linear) and therefore the contribution of each ion will add-up. The diffusion constants for different 
ions are of the same order of magnitude (for Na"*", K"*", CI", Ca^"*"), so no particular ion would be 
expected to dominate. 
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Tables 



n 


A6„ 


Tn is) 




an 


1 


4.50x10^ 


7.96x10- 


12 


0.10 


2 


4.00 X 10^ 


15.92X 10 


Q 


0.15 


3 


2.0 xlO^ 


106.1x10 


-6 


0.22 


4 


4.5x10^ 


5.305x10 


-3 


0.00 



Table 1: Parameter values for the parametric model of Gabriel et al. (fTS]) (see Eq.O. 



Figure Legends 



Figure 1: Illustration of the two main physical phenomena involved in the genesis of local field 
potentials. A given current source produces an electric field, which will tend to polarize the charged 
membranes around the source, as schematized on the top. The flow of ions across the membrane of 
the source will also involve ionic difi'usion to re-equilibrate the concentrations. This diffusion of ions 
will also be responsible for inducing currents in extracellular space. These two phenomena influence 
the frequency filtering and the genesis of LFP signals, as explored in this paper. 

Figure 2: Monopole and dipole arrangements of current sources. A. Scheme of the extracellular 
medium containing a quasi-dipole (gray) representing a pyramidal neuron, with soma and apical 
dendrite arranged vertically. B. Illustration of one of the monopoles of the dipole. The extracellular 
space is represented by cellular processes of various size (circles) embedded in a conductive fluid. 
The dashed lines represent equipotential surfaces. The ab line illustrates the fact that the extracellular 
fluid is linearly connex. 

Figure 3: Models of macroscopic extracellular conductivity compared to experimental measurements 
in cerebral cortex. The experimental data (labeled "G") show the real part of the conductivity mea- 
sured in cortical tissue by the experiments of Gabriel et al. (JT5). The curve labeled "E" represents 
the macroscopic conductivity calculated according to the effects of electric field in a non-reactive 
medium. The curve labeled "D" is the macroscopic conductivity due to ionic difi'usion in a non- 
reactive medium. The curve labeled "P" shows the macroscopic conductivity calculated from a reac- 
tive medium with electric-field efl'ects (polarization phenomena). The curve labeled "DP" shows the 
macroscopic conductivity in the full model, combining the eff"ects of electric polarization and ionic 
difi'usion. Every model was fit to the experimental data by using a least-square procedure, and the best 
fit is shown. The DP model's conductivity is given by ^ = i^o + Sji + ^ + with Kq = 10.84, 
Ki = -19.29, K2 = 180.35 and K3 = 52.56. The experimental data (G) is the parametric Cole-Cole 
model (fTO|) which was fit to the experimental measurements of Gabriel et al. |T3). This fit is in 
agreement with experimental measurements for frequencies larger than 10 Hz. No experimental mea- 
surements exist for frequencies lower than 10 Hz, and the difi'erent curves show diff'erent predictions 
from the phenomenological model of Cole-Cole and the present models. 



Figure 4: Simulation of 1/f frequency scaling of LFPs during wakefulness. A. LFP recording in the 
parietal cortex of an awake cat. B. Power spectral density (PSD) of the LFP in log scale, showing 
two different scaling regions with a slope of -1 and -3, respectively. C. Raster of 8 simultaneously- 
recorded neurons in the same experiment as in A. D. Synaptic current calculated by convolving the 
spike trains in C with exponentials (decay time constant of 10 ms). E. PSD calculated from the 
synaptic current, shown two scaling regions of slope and -2, respectively. F. PSD calculated using 
a model including ionic diffusion (see text for details). The scaling regions are of slope -1 and -3, 
respectively, as in the experiments in B. Experimental data taken from ref. (fTTI) : see also ref. (4) for 
details of the analysis in B-D. 



Figure 5: Simulation of more complex frequency scaling of LFPs during slow-wave sleep. A. Sim- 
ilar LFP recording as in Fig. |4K (same experiment), but during slow-wave sleep. B. Raster of 8 
simultaneously-recorded neurons in the same experiment as in A. The vertical gray lines indicate 
concerted pauses of firing which presumably occur during the "down" states. C. Synaptic current 
calculated by convolving the spike trains in B with exponentials (decay time constant of 10 ms). D. 
Power spectral density (PSD) of the LFP in log scale, showing the same scaling regions with a slope 
of -3 at high frequencies as in wakefulness (the PSD in wake is shown in gray in the background). 
At low frequencies, the scaling was close to 1 //" (gray line; the dotted line shows the 1 // scaling 
of wakefulness). E. PSD calculated from the synaptic current in C, using a model including ionic 
diffusion. This PSD reproduces the scaling regions of slope -2 and -3, respectively (gray lines). The 
low-frequency region, which was scaling as 1 // in wakefulness (dotted lines), had a slope close to -2. 
Experimental data taken from ref. (fTTI) . 
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